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Abstract 

We present a condition which guarantees spatial uniformity for the asymptotic behavior of the solu- 
tions of a reaction-diffusion PDE with Neumann boundary conditions. This condition makes use of the 
Jacobian matrix of the reaction terms and the second Neumann eigenvalue of the Laplacian operator on 
the given spatial domain, and replaces the global Lipschitz assumptions commonly used in the literature 
with a less restrictive Lyapunov inequality. We then present numerical procedures for the verification 
of this Lyapunov inequality and illustrate them on models of several biochemical reaction networks. 
Finally, we derive an analog of this PDE result for the synchronization of a network of identical ODE 
models coupled by diffusion terms. 

1 Introduction 

Spatially distributed system models are essential for many fields of science and engineering. In cell biology, 
gradients of protein activities organize signaling around cellular- structures and provide positional cues for 
important processes, such as cell division [1]. One of the theories for spatial organization and pattern for- 
mation is based on diffusion-driven instability [2, 3], which has been a subject of intense study as surveyed 
in [4-6]. This phenomenon occurs when one of the higher spatial modes in the reaction-diffusion partial dif- 
ferential equation (PDE) is destabilized by diffusion, thus causing nonuniformities to grow. Understanding 
when the solutions of a reaction-diffusion PDE exhibit uniform behavior is an important problem, because it 
rules out diffusion-driven instabilities and justifies a simpler ordinary differential equation (ODE) modeling. 
The standard approach to proving spatial uniformity in the literature is to establish exponential decay of ini- 
tial nonuniformities by using global Lipschitz bounds on the vector field representing reaction terms [7-10]. 
In the first part of this paper, we study the reaction-diffusion PDE: 

^=f{x) + DV^x, (1) 

subject to Neumann boundary conditions and other technical assumptions detailed in Section |2J and give 
a condition for uniform behavior of the solutions that does not rely on a global Lipschitz assumption on 
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/(x). Instead, our main result (Theorem [T) requires that a Lyapunov inequality be satisfied by the matrix 
J{x) - A2D, where 

J(.x) := — — (2) 
ox 

is the Jacobian and A2 is the second Neumann eigenvalue of the operator L = -V^x on the given spatial do- 
main. Even when the global Lipschitz condition of [7-10] holds, our result can achieve orders of magnitude 
improvements over the estimates obtained from this Lipschitz bound (see Example |2]for a comparison). 

In the second part of the paper (Section O, we parameterize 7(x) with constant matrices and develop 
procedures to verify the Lyapunov inequality employed in Theorem [T] The first procedure, described in 
TheoremEl incorporates /(x) within convex and conic hulls of constant matrices and derives a linear matrix 
inequality (LMI) [11] for the vertices. The second procedure, presented in Theorem [3j studies a special 
convex set and reduces the dimension of the LMI in Theorem |2l For reaction networks that exhibit special 
structures, the LMI in Theorem[3]is also amenable to analytical feasibility tests. One such test is illustrated 
in Example [2]on a variant [12] of Goodwin's model [13] for oscillations in enzyme synthesis. In Example 
|3j we study a model by Goldbeter [14] for circadian rhythms and investigate the feasibihty of the LMI 
numerically. 

In a recent study [15], we gave conditions for the stability of the spatially uniform fixed point for 
reaction-diffusion systems where the reaction terms exhibit a cyclic structure. In the present paper we 
do not restrict ourselves to cyclic reactions and, more importantly, we do not require that the attractor be a 
fixed point. Indeed, the reactions in Examples |2] and [3] exhibit limit cycles and Theorem [J guarantees spatial 
uniformity of the oscillations rather than stability of a fixed point. 

In the third part of the paper (Section |4ll, we derive an analog of Theorem [J for a finite number of 
identical ODEs coupled via diffusion-like terms [16]. For ODEs, the equivalent of spatially uniform behavior 
is synchronization, on which a large literature exists as reviewed in [17]. Our main result (Theorem |4ll in 
this part employs the same condition as Theorem [TJ where A2 now represents the second smallest eigenvalue 
of the Laplacian matrix for the graph describing the coupling of the subsystems. The proof of this result 
exploits properties of the Laplacian matrix that aie analogous to those of the Laplacian operator employed 
in Theorem [T] In Example |4l we make a connection between Theorem |4] and the incremental passivity 
approach to synchronization employed in [18]. 



2 Uniform Behavior in Reaction-Diffusion PDEs 

Let Q be a bounded domain in W with smooth boundary dQ., and consider ([T]) where x € R", /(•) is a 
continuously differentiable vector field, and V^x := [V^xi •••V^x,,]^ is the vector Laplacian. In a typical 
reaction-diffusion system, D € M."^" is a diagonal matrix of diffusion coefficients dj for each species / = 
1, • • • ,n; however, in the derivations below, we take D to be an arbitrary real matrix for further generality. 
We assume Neumann boundary conditions: 

Vx,(^)-n(^) = y^edQ, i=\,---,n (3) 

where ^ represents the spatial variable and h is a vector normal to the boundary dQ.. Well-posedness of ([B- 
(O is not emphasized in this paper; the reader may refer to [19, Chapter 7.3] for conditions that guarantee 
existence of classical solutions to reaction-diffusion PDEs. 
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To establish a condition under which solutions x{t,^) exhibit uniform behavior over the spatial domain 
Q, we denote by: 

7t{v} := V- V (4) 
the deviation of a function v = v(^) from its average: 

In the derivations below, we also use the ^2(0) inner product: 

{u,v}Mn):= f u\m^)d^ (6) 
Jn 

and norm: 

l|v||L2(n) yj{v,v)L2ia)- (7) 

We let = /ti < A2 < •■■ < < ■•■ denote the eigenvalues of the operator L = -V^ on Q. with Neumann 
boundary condition: 

Lcpt{€) = AkMa VM^)-m = o^^eda, (8) 

and make use of the second smallest eigenvalue, A2, in our main result: 

Theorem 1. Consider the reaction-diffusion system (CJ-fEl) and let A2 be the second smallest eigenvalue of 
the operator L - -V^ on Q. with Neumann boundary condition as in (H]). If there exists a convex set X c R", 
a matrix P = > 0, and a constant e > such that 

P{J{x)-A2D) + {J{x)-A2Df P <-eI VxeX (9) 
PD + D^P>0, (10) 

then, for every classical solution x(t,^) : [0, 00) x Q — > ?(, 

Mx{t,m\Mn)^0 (11) 

exponentially as t ^ 00. □ 

The second Neumann eigenvalue A2 is a measure of the well-connectedness of the spatial domain. In- 
deed, of all sets of given volume, A2 is maximized by the ball [20]. In situations where A2 is not easily 
calculable for the given domain Q, Theorem [U can be applied with a lower bound on A2 at the cost of mak- 
ing (l9l)-(fT0l) more restrictive. A commonly used lower bound on A2 was derived for the Laplacian operator 
by Cheeger [21], and extended in [22] to Laplacian matrices of graphs. 

Othmer [7], followed by other papers [8-10], studied the reaction-diffusion system ([T])-© with D = 
diag{(ii, • • • ,d„}, and proved uniform behavior of the solutions under the condition: 

sup||7(.x)|| < A2min{di}. (12) 

Note that ([T2l ) implies ^ with P = I, which means that Theorem [T] incorporates Othmer's condition ([T2l ) as 
a special case. Assumption ^ of Theorem [T] is far less restrictive than ([T2l ). and is applicable to numerous 
practically important systems which do not satisfy global Lipschitz bounds. 
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Example 1. As an illustration of Theorem [T] consider the Fitzhugh-Nagumo model of neuron excitation 
and oscillations (see e.g. [23]), augmented here with diffusion terms: 



~dt 



-TT" = C\x\-\x\ + X2\ + d\V^Xi 



9X2 1 2 

--— = -{-Xi-bx2 + a) + d2y X2, c,b,d\,d2>0. 
ot c 



The Jacobian matrix: 



J{x) 



c{\-x\) c 



does not satisfy a norm bound; however, with A2d\ > c, conditions (l9l)-(fT0l) hold with: 

c 



P = 



(13) 
(14) 

(15) 



(16) 



To see the implications of Theorem [J for a linear reaction-diffusion system, we recall that the solutions 
of ([T]) with f{x) = Ax can be expanded as: 



x(t,^) = Y,(^kit)MO (17) 

k=l 

where, due to the orthogonaUty of the eigenfunctions (pk{^) in ([D, crtit) e R" obey the decoupled ODEs: 

cfk^(A-A,,D)crk. (18) 

Because the eigenfunction (f>i corresponding to /li = is constant, the k = I term in ([17] ) constitutes the 
average x, governed by x = Ax, and the decay of the remaining terms is guaranteed when the matrices 
A - AjcD, ^ = 2, 3, • • • , are Hurwitz. Conditions (l9l)-(fT0l) with J{x) = A in Theorem [T] imply the existence of a 
common Lyapunov function for these matrices, thus ensuring that they are indeed Hurwitz. 



Proof of Theorem [TJ We denote 

X :- 7r{x] 

where n{-} is as defined in (011, and note that x satisfies: 

^=n{fix)} + DV^x, 



(19) 



(20) 



where we have substituted 7r{ V^x,} = V^x,- because V^x,- d^ - j^^ Vx,- -hdS =0 from the Divergence The- 
orem and the boundary condition (O. We then select the functional: 



Vix) = -{x,Px)l,jo.), 



(21) 



where P is as in (l9l)-(fT0l). and obtain: 



V < {x,Pn{f{x)])L,(ci) + {x,PDV^x)L,(m- 



(22) 
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We note from (ITOl) that there exists a matrix Q such that Q^Q = \{PD + D^P). This means that: 
where y :- Qx. Integrating both sides of the identity 



y-(yiyyi) = \yyif+yi'^bi (24) 

over Q and noting that the left-hand side vanishes due to the Divergence Theorem and the boundary condi- 
tion ([3l), we obtain: 

r yiV^yid^^- [ IVytfd^. (25) 

Moreover, because J^yd^ - Qj^xd^ = 0, it follows from the the Poincare Inequality [20, Equation (1.37)] 
that: 

r \yyifd^>A2 f yjd^ (26) 
Jn Jn 

and, thus, (1251 ) and (l26l ) imply: 

{yi,^\)L2(a)<-Myi\\l,^nr (27) 
We substitute the inequality (|27]l in (|23]l, substitute back y = Qx and ^ ^(PD D^P), and obtain: 

{x,PDV^x)L2(a) ^ {y,^b)Mm < -My,y}L2(a) - -^2<x,PD^)L,(n). (28) 
Substitution of (1281 ) in (l22l) then gives: 

V < <^,P7r{/(^)})t,(n) - ^2<^,P£>^)L2(n)- (29) 
Next, we rewrite the first term on the right-hand side of (l29l ) as: 

{x,P7T{f{x)])L2ia) - J^^^p|/(x)-il£/(x)J^J^/^ (30) 

= jj^P{f{x)-f{md^ + jj''p{f{x)--^^Jj{x)d^y^ (31) 

- r fP{f{x)-f{x))d^ (32) 
Jn 

where, to obtain (|3TI ). we added and subtracted f{x) in (l30l ). To obtain (l32l ). we noted that the second integral 
in (|3TI ) is zero because the factor 

does not depend on ^, and because j^xd^ = 0. Substitution of (l32l ) in (l29l) then results in: 

V < <x, P(/(x) - fix) - A2Dx))L2(a) ■ (34) 
Finally, we use the Mean- Value Theorem [24] and write: 

fix)- fix) = 1 Jix + six-x))ix-x)ds. (35) 
Jo 

Substituting in (l34l ) and using Q, we obtain: 

V< r r fPiJix + sx)-A2D)xd^ds< f { -^fxd^ds< ^ — V. (36) 

Jo Jn Jo Jn 2 /Imax(P) 

Inequality (l36l) proves exponential decay of the functional V(x) defined in (|2TI ). from which the conclusion 
of the theorem follows. □ 



3 Constant Matrix Parameterizations of the Jacobian 



We now present a procedure to verify ^ by bounding the Jacobian J(x) within a set which is parameterized 
by constant matrices. Examples of such parameterizations include the convex hull: 



conv\ 



4Zi,---,Zg} = [eiZi + --- + egZg\9i + --- + 9^ = l, 9i>0 i=l,--- ,q], (37) 
and the conic hull: 

cone{S\,--- ,Sm) = {coiS i H \-a>,„S,n \ coi>0 i = I,--- ,m}. (38) 

When J{x) belongs to the sum of these two sets, © can be replaced with the constant matrix inequalities 
(IlOll-diB below: 



Theorem 2. If there exist constant matrices Z\,--- ,Zq and 5 1 , • • • , 5 ,„ such that 

J{x) e conv{Zi,- ■ ■ ,Zg] + cone{S I,-- ■ ,S m] VxeX, (39) 

then a matrix P = satisfying: 

P{Zk-A2D) + {Zk-A2DfP<0, k^\,---,q (40) 
PSk + SlP<0, k=\,---,m (41) 

also satisfies &for some 6 > 0. If the image of X under /(•) is surjective onto 

conv{Z\,--- ,Zq] + cone{S \,--- ,5,,,}, 

then the converse is also true; that is, (|9l) with e> implies l\40i-l\41i. □ 

The proof is routine and is given in the Appendix. Theorem |2] is useful because the inequalities (ITOl ). 
(RTO and (BTl ) are linear in the variables P = P^ >0 and e > and, thus, the conditions of Theorem [1] can be 
checked with efficient numerical tools available for linear matrix inequalities [11]. Analytical conditions for 
the existence of common quadratic Lyapunov functions are also available for several classes of matrices [25]. 

In various examples of reaction networks, J{x) belongs to a convex set of the form: 

box{Ao,Ai,- ■ ■ ,A{} ^ {Ao + yiAi + ■ ■ ■ + y(A( \ <yi < I, /= 1, •••,€), (42) 

where A\, - ■ ■ ,Ac are rank-one matrices. Although Theorem |2] is applicable to the matrices Z\,--- ,Zq corre- 
sponding to the vertices of the set (l42l ). this application involves q = 2'' vertices and may become intractable 
for large £. Theorem |3] below, proven in the Appendix, gives an alternative test that uses only the matrices 
Aq, • • • , for verifying ®: 

Theorem 3. Suppose 

J{x)ebox{Aa,Au---,Ac}, (43) 
where A\,--- ,Ai are rank-one matrices and, thus, can be decomposed as: 

Ai = BiCj i=l,---,i (44) 
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with appropriately selected column vectors Bi, Ci e R". If there exists a matrix P = 'P^ > of the form: 



Pe. 



\ cjiER, i^ l,---J, 



(45) 



such that: 



A0-A2D 



B 



A0-A2D 



B 



P<0, 



(46) 



where B := [Bi ■ ■ ■ B{] and C := [Ci • • • Cf], then the upper-left block P = P^ >0 satisfies ^for some e>0. 
If, in addition, € = \ and the image of X under /(•) is surjective onto box{A(),Ai}, then the converse is also 
true; that is, ij® holds with a matrix P = P^ > and a constant e > 0, then there exists qi > such that 



P = P^ >Qin^ satisfies g^. 

Remark 1. In applications, it may be preferable to search for a fully diagonal matrix P satisfying 
instead of a block-diagonal P as in (|45]) . Although this restriction may add conservatism, it has the following 
advantages: 

i) Condition dTOl ) in Theorem [T] is satisfied for all diagonal and nonnegative D, and need not be checked 
separately when D has this form. 

ii) If the set box{Ao,A\, - ■ ■ ,Af} in Theorem[3]is augmented with cone{S i,--- ,S,„}, then the upper left 
n X n component P of the matrix P in (|45] ) must satisfy (BTl ) in addition to (|46l ). However, in the special case 
where S^-, ^ = 1, • • • ,m, are nonpositive diagonal matrices, (1411 holds for every diagonal P > and, thus, it 
is sufficient to check (l46l ) with a diagonal P>0. 

m) Likewise, if some of the matrices A, in (|43] ) are diagonal and nonpositive, the corresponding columns 
Bi and C, can be omitted in constructing the matrix 



A0-A2D B 
-I 



(47) 



thus reducing the dimension of the problem (|46l ). 

iv) In several practically important examples, analytical tests ai^e applicable to check the existence of a 
diagonal solution to the Lyapunov inequality (l46l ). Matrices ^ for which a diagonal P>0 satisfying 

(48) 



P^ + J{^P<0 



exists are termed diagonally stable [26], and have been fully characterized in dimension three [27] and di- 
mension four [28]. For higher dimensional matrices, diagonal stability tests have been derived by exploiting 
special sparse structures, such as a cyclic structure and its variants studied in [29, 30]. If the matrix Wh 
conforms to one of these structures, the existence of a diagonal solution P to (|46l ) can be checked with 
simple algebraic conditions. Analytical conditions ai^e indeed important in applications, because they reveal 
which system properties and which parameters are critical for Theorem[T]to hold. □ 
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Example 2. In [31], Jacob and Monod gave a molecular description of how certain metabolites regulate 
their production by repressing enzymes necessary for their synthesis. Following this description, Goodwin 
[13] proposed a differential equation model and studied its oscillatory behavior. A variant of Goodwin's 
model [12], augmented here with diffusion terms, is: 



dxi 
'df 

dx2 
dt 



Ki +X3 
-02X2 + biX\ + d2^^X2 

1^3X3 2 

+ Z72X2 +(i3V X3, 



(49) 



^3 + -^3 

where all parameters are positive and x\, X2, X3 denote, respectively, the concentrations of the messenger 
RNA, enzyme and product. 

To inspect condition ^ of Theorem [T] on the set X - 



l^Q, we study the Jacobian matrix: 



-a\ 
bi -a2 
b2 



-b3{X3) 



-«3te) 



and note that 



V3K3 



b^ix-i) :-- 



(^3+^3)2 (^1+^3)2 

lie in the bounded intervals [0, ■^] and [0, -p], respectively. This means that 

J{x) € box[AQ,Ai,A2}, 



(50) 



(51) 



(52) 



where 





-a\ 
















Vl 1 















bi 


-a2 





Ai = 











A2 = 
















b2 
























Vi 



We decompose Ai as Ai - B\C^ with Bi - [--p- 0]^ and Ci = [0 1]^, and construct the matrix: 



A0-A2D Bi 
C\ -I 



-{ax+A2dy) 

b\ -{a2 + A2d2) 



--^ 








b2 







-A2d3 
1 






-1 



(53) 



(54) 



where we have omitted A2 in view of item (iii) in Remark[T] Because the matrix (l54l ) has a cyclic form, the 
secant criterion derived in [29] is applicable, and states that diagonal stability of (l54l ) is equivalent to the 
condition: 



bib2Vi 



K7{a\ + A2d\){a2 + A2d2)A2d3 



< sec(7r/4)^ - 4. 



(55) 



We thus conclude from Theorem [3] and Remark [T] that, if the parameters of the model ( |49l ) are such that ( |55l ) 
holds with A2 calculated from the domain O, then Theorem [T] guarantees spatial uniformity of the solutions. 
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Note that, in this example, is bounded and, hence, condition ([T2l) of [7] is applicable. With the 

following set of parameters from [12]: 



a^=a2 = bi=b2 = 0.01, Vi = 9,V3 = Ki=K^ = 1, 



(56) 



([121 ) stipulates: 



A2mmdi > sup ||/(x)|| - 9.0554, 



(57) 



where the sup is achieved when = 0. To compare this condition to (1551) . we note that, for the same 
parameter values, 

bib2Vi 9-10-4 



Kj{ai + A2di)(a2 + A2d2)A2d3 (0.01 + /l2min,Ji)^'^2min,J,' 
which implies that (|55] | holds if the upper-bound in ( [58] ) is less than 4; that is, if: 

/l2min<i,- > 0.05435. 



(58) 



(59) 



The estimate (1591 ) is obtained using the upper-bound ( [58]) . which is achieved only when the diffusion co- 
efficients are identical. For nonidentical diffusion coefficients, condition ([55] ) leads to even more dramatic 
improvements over the conservative estimate (|57] ). In this example, a direct application of Theorem|2]to the 
vertices of the set ( [52] ). without insisting on a diagonal solution P, gave an insignificant improvement over 
(|59] ): With d\ - d2 = dj, - d, we numerically obtained the bound A2d > 0.05425. 

Example 3. As a further illustration of Theorem [3] we consider a model of Drosphila circadian rhythms, 
proposed in [14] and further studied in [32] for its dynamical behavior. When augmented with diffusion 
terms, this model is of the form: 



— - = k,M 



dM 
'df 
dP_Q 

dt 
5Pi 

dt 
dP2 

dt 

dt 



v„jM 
k,„ + M 



ViPo 



+ ■ 



+ dM'^ M 
V2P1 



Ki+Po K2 + P1 
ViPo V2P1 V3P1 



V4P2 



Ki+Pq K2 + P1 K3+P1 K4 + P2 

V3P1 V4P2 , „ . , V,,P2 



+ dp^V^Pi 



K3 + P1 K4 + P2 

kiP2-k2PN + dp^V^PN, 



-k^Pi+kiPN- 



+ dp,V^P2 



(60) 

(61) 
(62) 
(63) 
(64) 



where Pi represents the concentration of the PER protein, with the indices / = 0, 1 , 2 denoting the degree of 
phosphorylation. Likewise, P^ is the concentration of PER in the nucleus and M is the concentration of the 
messenger RNA. 



We obtain the Jacobian matrix 











-MPo) 
(piiPo) 









MPi) 
-MPi)-<f>3iPi) 
hiPi) 




-(f>5(PN) 



MP2) 

-ki-MP2)-MP2) k2 

k\ -k2 



(65) 
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where 
Ux) : 



KiVi 



/ = 1,2,3,4, 05W:= 



Vmkn 



and note that, for ;c > 0, these functions lie in the bounded intervals [0,^,], with: 



4>j(x) :-- 



(kd + xf ' 



^i = ^ /= 1,2,3,4, 



Thus, for all {M,Po,Pi,P2,Pn) e : 



= 7—, (Pl = 7-, 95 = 



4nK', 



if « - 1 



l^Q, 7€&ox{Ao,Ai,--- ,^7), where 



Ao 






















































V, 
K, 































Vi 
Kl 






















kl 




























kl 


-k2 _ 




















A7 = 
















-I 













A,- 

















































-4> 

























































V3 
K3 








A4- 











V4 





A5 = 























Ki 



















V4 











































































































































A6- 

















Ai = 





























































(66) 



(67) 



(68) 



(69) 



(70) 



Using the following parameter values from [14]: 

n = A, V, = 0.76, K[ - 1, k, = 0.38, ki - 1.9, k2 = 1.3, Vi = 3.2, V2 = 1.58, ^3 = 5, ^ = 2.5, 
=2,K2 = 2, = 2,K^ = 2, v,„ = 0.65, Vd = 0.95, k^ - 0.2, k,„ - 0.5, (71) 

and assuming identical diffusion coefficients, denoted by d, we applied the procedure outlined in Theorem 
[3] and numerically determined bounds for A2d using the MATLAB software CVX [33]. The linear matrix 
inequality (|45])-(|461) was feasible with a fully populated matrix P when A2d > 0.4590, and with a diagonal 
matrix P when A2d > 0.5393. 

Note that, in our parameterization, we took advantage of the repetition of the nonlinearities in 0i , • • • ^4 in 
(|65] ). and employed the matrices A 1, • • • ,A4, each representing two occurrences of the same nonlineaiity. The 
alternative approach of overparameterizing with one matrix for each occurrence would lead to conservative 
results. Indeed, a repetition of the numerical experiment described above with Ai,--- ,A4 split into two 
matrices each, gave the conservative feasibility region A2d > 1 .7892 with a fully populated P (compare to 
A2d > 0.4590 above), and A2d > 1.7943 with a diagonal P (compare to A2d > 0.5393 above). 
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4 Synchronization in a Network of ODEs with Diffusion-Like Coupling 



We now derive an analogous result for a network of identical ODE models that are interconnected according 
to an undirected graph: 

+ D 2] k = \,---,N, (72) 

where x'^ eR", NkQ{l,--- ,N] denotes the set of nodes adjacent to node k in the graph, and D is allowed to 
be an arbitrary real matrix as in Section |2] Denoting by X the concatenated vector: 

X^lx'^.-.x^'^f, (73) 

and by L = (/,• y) e R^^^ the graph Laplacian matrix [34] : 

kj = 

we rewrite ( [72b in the compact form: 



number of nodes adjacent to node / if / = j 

-I if ii= j and j e Ni (74) 
otherwise, 



X = F(X)-(L®D)X, (75) 
where "®" represents the Kronecker product, and 

FiX):=[f{x'f---f{x''ff. (76) 

We let = /li < A2 < ■•■ < Af^ denote the eigenvalues of the Laplacian matrix, and show that the compo- 
nents x'^(t) in ( [72] ) synchronize if A2 is such that (|9l)-([T0l) hold as in Theorem [T] 

Theorem 4. Consider the interconnected system d75D-d76D, and suppose f|9l)-diOD hold with a matrix P = 
> and a constant e > on a convex set X c R.". Then, every forward-complete solution X{t) = 
[x^{t)^ ■■■ x^{t)^] ^ that remains in has the property that, for any pair {k,j)^{\,--- ,N\x{\,--- , N), 

x\t)-xj{t)^0 (77) 

exponentially as t —> 00. □ 

Proof of Theorem m In this proof, we make repeated use of the property: 

(A(8)S)(C(8)D)-(AC)®(BD), (78) 

which holds whenever the matrices are of compatible dimensions to form the indicated products. We also 
recall that the Laplacian matrix L satisfies: 

LIn - (79) 

where 1/v denotes the Nxl vector of ones. Since 1/v is an eigenvector that corresponds to the eigenvalue 
A\ = 0, it follows that 

y'^Ly>A2y'^y Vj±1a,. (80) 
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Likewise, denoting by /„ the nxn identity matrix, we get the inequaUty: 

y^(L®In)y>A2y^y Vjil^®/,,, (81) 

which is the discrete analog of the Poincare Inequality (l26l ) used in the proof of Theorem [T] 
Mimicking ^ and (fT9l) . we define: 

x--^{x^ + --- + x^) = ^{ll^In)X, X:=In^x, (82) 



and 



x''--x^-x, X-X-X. (83) 



It follows from this definition that = and, thus, for any matrix M with n rows. 



N 



X^(1a,(8)M) = ^x*^^M = 0. (84) 
k=l 

The dynamics of X are given by: 

it = F{X)-X-{L^D)X 

= F{X)-X-{L®D)X, (85) 

where the second equation follows by substituting X - X + \i^®x and by noting from ( |78] ) and ( 1791 ) that 
(L(8)D)(1a,®^)= (LlAr)®(D^) = 0. 

We introduce the Lyapunov function V - ^X^ {In®P)X and note that it satisfies: 

y = X^{1m®P){F{X)-X)-X^{In®P){L®D)X 

= X^{1n®P){F{X)-X)-X^{L®{PD))X. (86) 

Because L is symmetric, the following identity holds: 

(L ® (PD)) + (L ® (PD))^ = L®{PD + D^P). (87) 

As in the proof of Theorem[T] we define Q such that Q = \{PD + D^P) and obtain: 

Z^(L ® (PZ)))Z = Z^(7^ ® e^)(L ® /„)(/^ ® Q)X - y'^{L®In)y, (88) 

where 

y:-(/;v®G)^. (89) 
It then follows from ^ and ([H) with M = that 

/(lw®4) = X^(/iv® e^)(lw®4) - ^""(liv® G^) - 0, (90) 

which means y 1. In® In and, thus, the inequality (ISTl) above is applicable. Using (ISTl) . (1881 ) and ( [89l ). we 
obtain: 

X^{L®{PD))X = y^{L®In)y > A2y^y = A2X^ (Jn®(PD))X = A2Y^y^^PDx!' . (91) 
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Substituting this inequality back in (1861) . we get: 

N 

V < X^(In ® P)(F{X) - X) - ^2 2 x^'^PDx^. (92) 

k=\ 

We next add and subtract F{X) - \iq®f{x), and rewrite ( |92l ) as: 

N 

k=i 

N 

= it^ilN ® P){FiX) - FiX)) + X^ilN® Pifix) - ^)) x'^^PDx'' 

k=i 

N 

= X^ilN®P)iF(X)-FiX))-A2Y,^''^PDx'', (93) 

k=i 

where the second equation follows from (1781 ) and the third equation follows from (l84l) with M = P{f{x) - x). 
Expanding the first term in ( |93] ) as a summation, we obtain: 

N N 

V<J] ^^P{f{^) - m) x'^'PDx^. (94) 

k=i k=i 

Finally, an application of the Mean- Value Theorem (1351 ) yields: 

N M 

V<y\ x'''P{J{x + sxi)-A2D)x'ds < -^X^X < -^-^V, (95) 

which concludes the proof. □ 

Remark 2. In Theorem |4l we assumed an undirected graph to give an exact analog to the reaction-diffusion 
PDE result of Theorem [T] However, with the additional condition that the product PD be symmetric, it is 
not difficult to extend Theorem |4] to a directed graph, where L is not symmetric and is restricted only by 
LIn = 0. In this extension, ( [87] ) must be replaced with: 

(L ® (PD)) + {L® {PD)f = {L + L^)® (PD) (96) 

which holds because PD is symmetric, and A2 must be redefined as the largest number such that dSOl ) holds. 
This definition of A2 was introduced in [35] as the "algebraic connectivity" of a directed graph, and employed 
in [36] to obtain a synchronization result over directed and weighted graphs. □ 

It is important to note that the Lyapunov inequalities ([9l)-([T0l) used in Theorems [1] and [4] imply a con- 
traction property [37] for the family of vector fields {/(x) - ADx, A> A2]. Contraction properties, in various 
forms, have been employed in [18,36,38-41] to derive synchronization conditions for networks. We now 
make a connection between Theorem [4] and one of the results in [18]: 

Example 4. Stan and Sepulchre [18] studied the ODE model^^: 

.%k _ A ,Jc , D,f,.k\ , D..k 



X 



+ B(p{y'') + Bu'' (97) 
/ - Cx\ (98) 



'We follow a slightly different notation than [18] for consistency with Theorem|4| 
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x'' eW, u'' eR, y'' eR, k = I,-- - ,N, coupled by the feedback law 



u = -Ly, (99) 

where u:-[u^ ■ --u^]^ and y :-[y^ ■■ ■y'^V, and pursued incremental passivity arguments to prove synchro- 
nization of the subsystems. 

We now show that Theorem 0] is applicable to (|97])-(l99l) when the following hypotheses, adaptec^l from 
those in [18], hold: 

HI. There exists a constant y such that (p'{y) <y, Sy e R. 

H2. The triplet {A + y*BC,B,C) is strictly positive real [42]; that is, there exits P = pT >Q such that: 

P{A + y*BC) + (,A + y*BCfP<0 (100) 
PB = C^. (101) 

H3. Llf^ = and the largest number, A2, such that ( ISUl ) holds, satisfies: 

A2>y-y*. (102) 

To apply TheoremU note that system is of the form with 

f{x)=Ax + B(f){Cx) D = BC. ( 1 03) 

From HI, we conclude that the Jacobian J(x) is as in ( [39l ). with 

Zi = A + yBC S 1 = -BC. (104) 

Noting from (llOll l that 



PD = PBC = C^C > (105) 



and using ( 11021) . we obtain: 



P(Zi - A2D) + (Zi - A2DfP < PiZi - (r - y*)D) + (Zi - (y - y*)DfP (106) 

Substituting Zi from ( 11041 ) in the right-hand side of ( 11061 ) and using dlOOK we conclude that condition ( |40l ) of 
Theorem|2]holds. Likewise, (|4T] ) holds because P5 1 = -PBC = -C^C < 0, and Theorem |2] verifies condition 
^ of Theorem |4] on /V = R". Finally, noting from (1105b that (fTOl) also holds, and that PD is symmetric as 
stipulated in Remark[21 we conclude dTTl ) for all forward-complet^ trajectories. □ 



^Unlike [18], in HI, we assume that the nonlinearity (p(-) is difFerentiable. In H2, we strengthen the passivity assumption of [18] 
to strict passivity so that dlOOI l holds with strict inequality and, thus, Theorem|4]is directly applicable. With a slight modification of 
Theorem|4l it is indeed possible to remove the strictness condition and, instead, to assume observability of the pair (A,C) as in [18]. 
Finally, in H3, we remove the "balanced graph" assumption (1^^ = 0) employed in [18]. 
[18] indeed argues boundedness for the solutions of (|97|)-l|99b. using ideas from [43]. 
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Appendix: Proofs for Theorem |2] and Theorem |3] 

Proof of Theorem IS From for every x e ?(, there exist parameters 9i,-- ■ ,9q,(±)\,--- , aj,n >0,9i-\ 1- 

dq = I, such that: 

J{X) = 9iZi + ■■■ + 9qZq + OJiS\+--- + OJniSm- (107) 

Because 9\-\ \-9q- 1, we write: 

J{x) - A2D = 9x{Zx - A2D) + ■■■ + 9q{Zq - A2D) + a;iS 1 + ■ • • + w^S^, (108) 

from which it follows that a matrix P satisfying (I^Ut-lPTl) also satisfies ©. To prove the converse, we note 

from the surjectivity assumption that, for any set of parameters 9i,--- ,9q,(jD\,--- ,a>,„ >0, 9i-\ \-9q = I, 

there exits x e <Y for which (11081 ) holds. To see that Q implies (l40l ). pick 9k = I, 9i = ii= k, and a*; = / = 
1, • • • ,m in (11081) . To see that ^ implies (|4TI ). assume, to the contrary, that ^ holds, but (|4TI ) fails for some 
k, which means that there exists ^ e R" such that 

^V5t + 5[P)^>0. (109) 

Then, pick Ui-Oii^k and note from ( 11081 ) that the left-hand side of ^ is equal to: 

q 

co^iPSk + SlP) + 2 9k[P{Zk - A2D) + {Zu - A2DfPl (110) 

k=i 

Because < 0^ < 1, choosing oj^ > large enough ensures that 

l^cok{PSk + S IP) + g 9k{P{Zk - A2D) + {Zk - A2DfP^ ^ > 0, (111) 

which contradicts □ 
Proof of Theorems We rewrite (l45l)-(|46l) as 

Ar^- A'.n\''p PR + rn 

<0, (112) 



P(Ao-A2D) + (Ao-A2DfP PB + CQ 



QC^ + B^P -2Q 
where Q = diag{q\ ,q{}, and make use of the following lemma, proven separately below: 
Lemma 1. If there exists a diagonal ixi matrix Q>0 satisfying ( lii2D . then 

P{Ao + riBiCf + ■ ■ • + jfBeCj - A2D) + (Aq + nSjCf + • • • + jcBiC^^ - A2DfP < Vy; e [0, 1]. (113) 

When t = \, the converse is also true; that is, if Ml 3^ holds for every ji e [0, 1], then ([112i holds for some 
constant Q> 0. 

To conclude the first statement of Theorem [3j we note that (11 131 ) implies ^ with e > 0. To prove the second 
statement, we note from the surjectivity assumption that ^ with e > implies (11131 ). Since £ = 1, we apply 
the converse statement in Lemma \T\ and conclude that (II 121 ) holds for some constant Q > 0; that is, (l46l ) 
holds with ^ = 1 and !P as in (|45l). □ 
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Proof of Lemma[TJ Defining F to be a diagonal matrix with entries y, e [0, 1], we rewrite (11131) as 



x^[P{Aq - AiD + BrC^) + (Ao - A2D + BTC^fP]x < Vx ?t 0. 
We then define the new variable 



and rewrite (11141) as: 



3; := rax, 

P{Ao-A2D) + {Ao-A2DfP PB 
B^P 



X 

y 



<o. 



Next, we note from ( 11151 ) with y,- € [0, 1] that yi and Cjx are constrained by: 

1 



yi{Cfx) = -y^>yl 
Ji 



which means that: 



eiC] 



-2eie] 



>0 /=1,-- 



(114) 



(115) 



(116) 



(117) 



(118) 



where e,- is the /th unit vector in R^. Thus, (11141 ) is equivalent to the statement that (11161 ) holds for all xi^O, 
yi^O, satisfying (fTTSl) i=l,--- 

We now invoke the S-procedure [11] which states that, for symmetric matrices Tq,T\, - ■ • ,T(, 



fPo^ < for all ^ 7t satisfying fPi^ > / - 1, • • • , ^ 
if there exist q\ >0,--- ,qc>0 such that 

To+qiTi + --- + qfTc<0. 



(119) 



(120) 



Because the matrices in (II 161 ) and (11181 ) play the roles of Tq and T,, / = I,-- - in the S-procedure, we 
conclude that (11161 ) holds for allxi^O,yi^ 0, satisfying (11181 ) if 



P(Ao-^2^) + (Ao-/l2£>)'^/' PB 
B^P 



ce 

QC^ -2Q 



<0 



(121) 



for some diagonal Q>0. Finally, we note that (|121l) requires 2 > because, if Q contains zero diagonal 
entries, then the matrix in (|121l) also contains zero diagonal entries and, thus, cannot be negative definite. 
This concludes the proof of the first statement of the lemma, because inequality (11211) is identical to (I112I ). 

To prove the converse statement, we recall that, when £ =1, the S-procedure also states that ( 11191 ) implies 
(11201 ) for some q\ > 0, provided there exists ^0 such that ^qT\^o > 0. When £ = I, Ti defined in (11181) has 
the form: 

C 
-2 



(122) 



which means that ^^T\^o > indeed holds with the choice ^0 = [^0 ^(^^xq)]^, where xq is such that C^xq + 
0. Because To is as defined in (I116I ). we conclude from the S-procedure that if (11131 ) holds for every 
71 e [0, 1], then (11211) holds for some constant 2 > 0. Recalling that (11211 ) cannot hold if 2 = and that 
(11211) is identical to (11121 ). we conclude that (11121) must be true for some Q > 0. □ 
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